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Abstract 

We study the Lyapunov exponent for electron and phonon excitations, 
in pure and random Fibonacci quasicrystal chains, using an exact real 
space renormalization group method, which allows the calculation of the 
Lyapunov exponent as a function of the energy. It is shown that the Lya- 
punov exponent on a pure Fibonacci chain has a self-similar structure, 
characterized by a scaling index that is independent of the energy for the 
electron excitations, "diagonal" or "off-diagonal" quasiperiodic, but is a 
function of the energy for the phonon excitations. This scaling behav- 
ior implies the vanishing of the Lyapunov exponent for the states on the 
spectrum, and hence the absence of localization on the Fibonacci chain, 
for the various excitations considered. It is also shown that disordered 
Fibonacci chains, with random tiling that introduces phason flips at cer- 
tain sites on the chain, exhibit the same Lyapunov exponent as the pure 
Fibonacci chain, and hence this type of disorder is irrelevant, either in the 
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case of electron or phonon excitations. 
PACS: 71.23.Ft, 71.23-k, 63.50.+X 

1 Introduction 

The experimental discovery of quasicrystals/ and also the building of artificial 
multilayer structures by molecular beam epitaxy,^ have considerably stimulated 
the theoretical study of quasiperiodic systems. Quasicrystals have a deter- 
ministic aperiodicity that characterizes them as intermediate structures between 
periodic crystals and disordered materials, therefore being expected to display 
new behavior. There has been in particular, great discussion on the nature 
of the energy spectrum and eigenstates of electron and phonon excitations on 
quasicrystals. It is questioned whether the spectrum is absolutely continuous, 
pointlike or singular continuous, or correspondingly, if the states are extended, 
localized or critical. 

The Fibonacci chain is the simplest quasicrystal, a one-dimensional sys- 
tem where the site or bond variables take one of the two values A and B, 
and are arranged in a Fibonacci sequence. The Fibonacci chain can be con- 
structed recursively by successive applications of a substitution rule, A AB 
and B A, oi alternatively, by successive applications of a concatenation 
rule, Si = Si-i ® Si-2, Si being the Fibonacci sequence at iteration i. The 
quasiperiodicity of the Fibonacci chain is characterized by the golden mean 
T = (l + VS) /2 , which gives the ratio of the number of A and B units. 
Tight-binding electron and phonon excitations have been studied on a Fibonacci 
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chain, using mainly transfer-matrix^"^^ and real space renormalization group 
techniques. It has been found that the energy spectrum for these excita- 
tions is a Cantor set with zero Lebesgue measure, this result having in addition 
been proven^^ for the case of electronic excitations on a site Fibonacci chain. 
The spectra of the periodic approximants to the Fibonacci chain exhibit self- 
similarity in the band structure, with a scaling index that for the electronic 
excitations is independent of the energy, while for the phonon excitations it 
is a function of the energy. ^ The integrated density of states for the various 
excitations presents rich scaling behavior, with indices varying from the edge 
to the center of the bands. ^°~^^'^^The characterization of the eigenstates on a 
Fibonacci chain is a more difficult task, and it has usually been restricted to 
a few special energies on the spectrum, for which the states are found to be 
self-similar or chaotic. More generally, it has been found evidence for the states 
being neither extended nor localized in the usual sense. ^^-^^ 

The localization properties of the states can be studied through the cal- 
culation of the Lyapunov exponent 7, which characterizes the evolution of a 
wavefunction along the chain.^''"^^ The Lyapunov exponent is zero for an ex- 
tended or critical state, but is positive for a localized state representing then 
the inverse of the localization length. Delyon and Petritis^° have proved that 
the Lyapunov exponent for a class of binary quasiperiodic tight-binding chains, 
vanishes on the spectrum, which rules out the presence of localized states. The 
Fibonacci sequence does not however belong to this class of chains, and the 
characterization of the states in that case remains under discussion. A study 
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on localization lenghts of tight-binding electrons on a pure Fibonacci chain has 
been presented by Capaz et al.,^^ that found no evidence for localization of the 
states. 

Real systems have in general some disorder. Random quasicrystals, in the 
sense of a random tiling, have been considered^^ to explain the properties of qua- 
sicrystalline alloys. It is well known that disorder has pronounced effects on the 
transport properties of periodic systems, specially in one-dimension where all 
the states turn to localized whatever the amount of disorder. A striking property 
of quasicrystals is that they exhibit extremely high resistivities, which decrease 
with the amount of defects. The effects of some types of disorder on the elec- 
tronic spectra and wavefunctions of Fibonacci chains have been consider ed.^^~^^ 

In this work we study the Lyapunov exponent for electron and phonon ex- 
citations in pure and random Fibonacci quasicrystal chains. We consider elec- 
trons in a tight-binding model, with " diagonal" -site and "off-diagonal" -bond 
Fibonacci ordering, and phonons on a lattice with bond Fibonacci ordering. 
The disorder introduced is random tiling imposed on the substitution or con- 
catenation rules for construction of the Fibonacci chains. We use a real space 
renormalization group method, which allows the calculation of a wavefunction 
along the chain, for any given energy, and therefore enables the determination 
of its Lyapunov exponent. This method provides a simple and very efficient way 
of numerically calculating the Lyapunov exponent as a function of the energy, 
for large Fibonacci chains. The method has great similarities, but also impor- 
tant differences, as will be discussed, with that used by Capaz et aJ.^^ The 
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method is based on decimation, which is here applied either to substitution^^ 
or concatenation,^^ and implemented in the presence of disorder. 

In order to calculate an eigenstate, one needs to specify an energy on the 
spectrum. Since the spectrum of a Fibonacci chain is a Cantor set with zero 
Lebesgue measure, the probability of numerically specifying an energy on the 
spectrum is essentially zero. Hence any chosen energy will almost certainly 
correspond to a gap, and the calculated Lyapunov exponents are then associated 
to gap states. It is shown that the Lyapunov exponent for the gap states of the 
various excitations has a fractal struture, and we study its scaling properties. 
Prom these properties we obtain information on the Lyapunov exponent for the 
states on the spectrum of the Fibonacci chain, and therefore on the localization 
properties of the excitations. We study the Lyapunov exponent for both, tight- 
binding electrons and phonons, remarking that the Goldstone symmetry present 
in the later and absent in the former, may lead to important differences in the 
scaling properties of the two systems. 

The outline of the paper is as follows. In Sec. II we describe the tight-binding 
electron and phonon systems that are studied, and present the renormalization 
method used to calculate the Lyapunov exponent. In Sec. Ill we present the 
Lyapunov exponent for the various excitations on a pure Fibonacci chain, study 
its scaling properties and discuss localization, and finally consider the effects of 
disorder on the Lyapunov exponent. In Sec. IV we present our conclusions. 
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2 Renormalization approach 

The dynamics of tight-binding electron and phonon excitations on a Fibonacci 
quasicrystal chain, can be described by the generic equation 

(e„-i;)*„ = F„_i*„_i + y„*„+i. (1) 

For the electrons, denotes the amplitude of the wavefunction at site n, 
corresponing to energy E, £„ is a site energy, and Vn is the hopping amplitude 
between site n and n + 1. For phonons, represents the displacement from 
the equilibrium position of the atom at site n, E = mu)^, to being the phonon 
frequency and m the atom mass, £„ = Vn-i + Vn, and Vn is the spring constant 
connecting sites n and n + 1. This latter model describes equally well spin 
waves on an Heisenberg ferromagnet at zero temperature, replacing the spring 
constants by exchange constants, and mo;^ by the spin wave frequency w. We 
note the Goldstone symmetry present in the phonon system, which imposes a 
correlation between the site e„ and the coupling Vn parameters, that is not 
present in the electron system. 

The various Fibonacci quasicrystal models are defined as follows. For elec- 
trons, the "diagonal" model is obtained from (1) by setting, = 1 and e„ = ea 
or eb, according to the Fibonacci sequence, and the "off-diagonal" model is 
obtained from (1) by setting, = and Vn = Va ov Vb, according to the 
Fibonacci sequence. The model for phonons is obtained from (1) with the cou- 
plings Vn = Va or Vb , arranged in the Fibonacci sequence. 
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The disordered Fibonacci chains are built by introducing random tihng in 
the substitution rule for construction, 

B ^ A, 

A AB, probability p, (2) 

A — > BA, probability 1 — p, 

in each iteration i, starting with B, the two possibilities corresponding respec- 
tively to direct and inverse substitution, or they are built by introducing random 
tiling in the concatenation rule for construction. 

Si = Si-i (g) Si-2, probability p, (3) 
Si = Si-2'SiSi-i, probability l—p, 

starting with So = B and Si = A, the two possibilities corresponding respec- 
tively to direct and inverse concatenation. Random tiling on substitution or 
concatenation generates, at each iteration, an identical set of disordered Fi- 
bonacci chains, though throug a different sequence of preceeding chains (e.g. 
A ^ AB ^ ABA ABABA, by substitution, vs, A ^ BA ^ AAB 
ABABA, by concatenation). 

The method that we use to calculate the Lyapunov exponent is based on 
the fact that the wavefunction at the Fibonacci sites n = n{i) = -Fj+i, 
given by Fj+i = Fj + F,_i with Fi = Fq = 1, can be easily obtained via a 
real-space renormalization group transformation, which consists in eliminating 
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appropriated sites on the chain, so that a chain similar to the original one is 
obtained, with a rescaled length and renormalized parameters. Under successive 
decimations one carries the system through larger length scales separating the 
sites. For the Fibonacci chain it is possible to deduce an exact renormaliza- 
tion transformation for the parameters £„ and Vn, the rescaling factor, under 
which the system is self-similar, being equal to r. After i iterations, the renor- 
malization transformation takes, for example, Va to V^'^ , which represents the 
renormalized interaction between two sites that are a distance r' apart, mea- 
sured in units of the original lattice spacing. The Fibonacci sites n{i) become 
the first neighbours of the end site n = 0, at each iteration i. Now, writing (1) 
as a recursion relation for the wavefunction, and fixing the "free-end" boundary 
condition V-i = 0, one gets, = [{eq — E) /Vb] \E'o- The wavefunction at 
the consecutive Fibonacci sites n{i), can therefore be obtained in terms of the 
parameters under successive renormalization iterations i, through 



fixing 4*0 (e.g. \l/o = 1). The Lyapunov exponent 7 is then calculated from the 
wavefunction (4), given that 



and Xn = r' for n = n{i). In the work of Capaz et al. the localization of the 
wavefunction \1/ is studied following the behavior of the coupling V under suc- 
cessive renormalizations, and not through the evolution of the wavefunction (4), 
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which also involves the parameter e. Although the behaviour of ^ is mainly de- 
termined by V, the complete expression should be used. Furthermore, in that 
work a small imaginary part rj is added to the energy E, which produces an 
artificial decay of the coupling V, that alters the actual evolution of the wave- 
function, and consequently interferes in the study of localization and evaluation 
of the Lyapunov exponent. 

Now we present the decimation techniques used to obtain the renormal- 
ization transformation of the parameters eo and Vo, for chains constructed by 
substitution or concatenation. 

A. Substitution chains 

The renormalization transformation of the parameters is obtained by elimi- 
nating sites in such a way as to reverse the substitution procedure in (2).^'* In 
order to build the transformation one needs to consider an expanded parameter 
space, for the various excitations, where the bonds Vn assume two different val- 
ues, Va and Vb, arranged in a Fibonacci sequence, and the site energies e„ may 
assume three different values, depending on the local environment of n, if 
Vn-\ = V„ = Va, 60 if Vn-\ = Va and K = Vb, £7 if K-i = Vb and Vn = Va- 
A choice of the initial parameters Va , Vb, Sa, £/3, £7, casts the problem into 
the model for electron excitations, "diagonal" {Va = Vb, Sa = ^-y ^p) OJ" 
"off-diagonal" {Va 7^ Vg, = = s-y), or phonon excitations {Va 7^ Vb, 
Ea = 2Va, £/3 = £7 = Va + Vb)- The reversal of rule (2) is achieved through 
the elimination of (3 — sites , corresponding to direct substitution, or 7-sites, 

9 



corresponding to inverse substitution. The resulting renormalization equations 



are: 

i) direct substitution, 



e'r' = e^'-V^''/{ef-E), (6) 
and for the end site n = 0, 



4+1) = - ^« V - i^) , v^J'+i) = y^i), (7) 

ii) inverse substitution, 

4^'^ = el:^-Vr/{e^;^-E), (8) 
4^+1) = s^-V^'/{e^^-E), 

and for the end site n = 0, 



-(i+l) 



V^^^ / (e'^) - e) , vt'^=V^;+'\ if ^0^'^=^, 



.(0 
-0 ' 



(i+1) 



if V^n 



(i) 



(9) 
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B. Concatenation chains 



The renormahzation transformation of the parameters is obtained by elim- 
inating the central site, after having performed concatenation, so as to reverse 
the concatenation procedure (3).^^ This leads to the following renormahzation 
equations, which are different for bond Fibonacci ordering, i.e. "off-diagonal" 
electrons and phonons, or site Fibonacci ordering, i.e. "diagonal" electrons. 

For the bond problem: 

i) direct concatenation. 



= et!^-yt''V{e^T'-E), (10) 
ii) inverse concatenation. 



= ^-Vri(fi^''-E), (11) 

with, = e^^'^^+ef, and the initial values, V^'^^ = Vb, V^^^^ = Va , = 

— ~ for "off-diagonal" electrons, and e'^^ = e^^ = Vq^^ = Vb, 
Sq^^ = Ei'^ = V^^'^ = Va, for phonons. 
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For the site problem: 
i) direct concatenation, 



^(i+l) 

^0 


= 4"- 






= 4-'.' 


- vt 









(12) 



E 



ii) direct concatenation, 



-(i+l) 
-0 

Ji+l) 



with the initial values, V^'^'> = T, e^^^ = e^^^ = sa - T'^/{eB-E), V^^^ 



(13) 



T"^/ {sb — E), and in i) e^'' = ea, ^i^^ = es, while in ii) e^'^ = £b, = sa- 

Considering the general case of a random Fibonacci chain, for a given prob- 
abihty of disorder p, we start with a specific disordered configuration, generated 
by (2) or (3), respectively for substitution or concatenation chains, and then 
iterate (6) —(9), (10) — (11) or (12) — (13), depending on the system studied, 
according to that configuration, in order to obtain the successive values for Vj}'' 
and £q\ This allows us to calculate the wavefunction provided by (4), at 
the successive Fibonacci sites, for a given energy E. For each probability p, we 
average the obtained wavefunction for E over many different disorder configu- 
rations. It is important to remark that when dealing with random chains, one 
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should first calculate the wavefunction for a specific disordered configuration 
and then average over configurations, instead of averaging the parameters over 
disorder at each step of the renormalization and then calculate the wavefunction 
with the averaged parameters. This latter procedure^^ will wash out important 
correlations in the system, and leads to different results depending on how the 
average is performed. The first procedure describes the physics more accurately. 

3 Lyapunov Exponent for Fibonacci Chains 

We now present the results concerning the Lyapunov exponent, calculated as 
a function of the energy, for the tight-binding electron, "diagonal" and "off- 
diagonal" , and phonon excitations on the pure and random Fibonacci chains. 
We consider first the case of pure chains, for which we study the scaling prop- 
erties of the Lyapunov exponent and their implications for the localization of 
states on the spectrum, and analyse afterwards the effects of disorder, of the 
kind of random tiling, on the Lyapunov exponent. 

As mentioned above, the wavefunctions that we numerically calculate cor- 
respond to gap states. Figure 1 shows the typical behavior of a wavefunction 
^n, at any chosen energy E, either for the electron or the phonon excitations 
on a pure Fibonacci chain. One observes that the wavefunction first oscillates 
over a certain length, and then grows exponentially. This behavior has mixed 
characteristics of an extended (oscillating) state and a localized (exponential) 
state. The length over which a wavefunction oscillates is a "memory" length, ^° 
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in the sense that beyond this length it loses memory of its initial phase. The 
exponential growth of the wavefunction is characterized by the Lyapunov expo- 
nent, which measures the inverse of a "locahzation length". We find that the 
"memory" length ^ and the Lyapunov exponent 7 are simply related, ^ w I/7. 
In figure 2 we present the Lyapunov exponent for the electron, "diagonal" and 
" off-diagonal" , and phonon excitations on the pure Fiboncci chain, calculated as 
a function of the energy. The exponent exhibits a rather nontrivial dependence 
on the energy, which has a clear correspondence with the associated density of 
states obtained by Ashraff and Stinchcombe^'*'^^ for the various cases, the finite 
values of the Lyapunov exponent corresponding to gap states, with the further 
a state is inside a gap the larger is its Lyapunov exponent. The Lyapunov ex- 
ponent exhibits a fractal structure, i.e. under dilation the same structure is 
revealed in a smaller scale, as can be seen by comparing the Lyapunov plots 
in figure 2 with those in figure 3. This structure is observed even in the very 
low energy range of the magnetic excitations, where 7 takes particularly small 
values, most probably due to the Goldstone symmetry. 

The scaling behavior of the Lyapunov exponent is studied through the vari- 
ation of the maximum exponent in a gap, 7max, versus the gap width, AEg.^^ 
We find that 

7„,ax ~ {^EgY, (14) 

where the scaling index 6, is independent of the energy for the electron exci- 
tations, "diagonal" and "off-diagonal", as shown in figure 4, but depends on 
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the energy for the phonon excitations, as figure 5 reveals, and it is shown in 
figure 6. We also find that the scaling index for the electron excitations de- 
pends on the quasicrystal site (e^, £b) or bond {Va,Vb) parameters, decreasing 
as the difference between the parameters increases, while the scaling index for 
the phonon excitations, varying with energy, also depends on the quasicrystal 
parameters {Va,Vb)- Our results for the electron excitations are in agreement 
with those obtained by Capaz et al.,^^ though their scaling indices differ from 
ours, probably due to the fact that they have calculated the Lyapunov exponent 
from the behavior of the coupling V alone and not from the evolution of the 
wavefunction ^, in (4), and moreover have introduced an imaginary part in the 
energy which influences the Lyapunov exponent, as discussed earlier. 

From the scaling expression (14), one obtains, for the various excitations, 
that 7max — > when AEg — > 0, implying that the Lyapunov exponent for 
wavefunctions on the spectrum, vanishes. We therefore have that the electron, 
" diagonal" or " nondiagonal" , and phonon excitations on a Fibonacci chain are 
nonlocalized. 

Let us now study the effects of disorder on the Lyapunov exponent. Disorder 
has drastic effects on the wavefunctions of one-dimensional periodic systems, 
localizing all the states. Figure 7 illustrates this fact, showing the Lyapunov 
exponent for phonon excitations on a random periodic chain, with couplings 
Va and Vb, as a function of the probability p of disorder, for various energies. 
One sees that the Lyapunov exponent increases with disorder, being also an 
increasing function of the energy. 
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For the random Fibonacci chains we considered disorder of the kind of ran- 
dom tiUng, introduced in the substitution or concatenation rule for construction 
of the chains. The resulting disordered chains differ from the pure Fibonacci 
chain in having a varying number of phason flips, located at certain points on 
that chain. By a phason flip it is meant a local rearrangement of tiles on the 
quasiperiodic structure, corresponding to a switch of the site, ea and £b , or the 
bond, Va and Vb, parameters. Using the cyclic property of the trace one can 
see that all those random tiling chains have the same spectrum, for the electron 
and the phonon excitations, as the pure Fibonacci chain. In the work of Lopez 
et aJ.,^^ on the effects of that kind of random tiling on the electronic excitations 
of a Fibonacci chain, it has however been found that the disorder affects the 
spectrum of the excitations. We think that this result is a consequence of the 
averaging of the parameters over disorder taken at each step of the renormal- 
ization in that work, which as already mentioned, loses important correlations 
in the system and introduces effects that depend on the averaging procedure 
used, corresponding in fact to different systems. On the other hand Naumis 
and Aragon,^^ considering electronic excitations, have also noted that phason 
flips located at certain points on the Fibonacci chain do not alter the spectrum 
of the excitations. 

We calculated the Lyapunov exponent for the electron and phonon excita- 
tions on Fibonacci chains with random tiling, as a function of the probability 
of disorder, for different values of energy. The results obtained are illustrated 
in flgure 8. We flnd that the disorder considered does not affect the Lyapunov 
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exponent either for the electron, "diagonal" or "off-diagonal", or the phonon 
excitations. The same result is obtained for disordered Fibonacci chains with 
random tiling either in the substitution or the concatenation rule for construc- 
tion of chains. The irrelevance of disorder found for the Lyapunov exponent of 
excitations on a Fibonacci chain is surprising, knowing the drastic effects that 
disorder has on the excitations on periodic chains. However, it should be noted 
that random tiling introduces a kind of bounded disorder, which has also corre- 
lations, and therefore might not be sufficient to produce localization of states. 
Furthermore, in contrast to the general case, it has been reported that there 
exist particular random potentials in one dimension that allow for extended 
states, those being described by an iterative procedure of construction.^'"'^^ Liu 
and Riklund^'* have found that other types of disorder, different from the one 
considered by us, produce localization of electronic excitations on a Fibonacci 
chain. 

4 Conclusions 

We have studied the Lyapunov exponent for tight-binding electron, "diagonal" 
and "off-diagonal", and phonon excitations in pure and random Fibonacci qua- 
sicrystal chains, using a real space renormalization group method. This method 
allows the calculation of a wavefunction along the chain, and the determination 
of the associated Lyapunov exponent, as a function of the energy, in a very 
efficient way for very long chains. We have found that the Lyapunov exponent 
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for the pure Fibonacci chain has a self-similar structure, being characterized by 
a scaling index that is independent of the energy for the electronic excitations, 
but depends on the energy for the phonon excitations. The scaling properties of 
the Lyapunov exponent, imply that it vanishes on the spectrum for the various 
excitations. We therefore have that the electronic and phonon excitations are 
nonlocalized on the Fibonacci chain . Considering random Fibonacci chains, 
we calculated the Lyapunov exponent as a function of the probability of dis- 
order, and found that the disorder introduced, of the kind of random tiling, 
does not affect the Lyapunov exponent, which takes the same value as for the 
pure Fibonacci chain whatever the degree of disorder, either for the electron 
or for the phonon excitations. The random tiling considered generates in fact 
chains that are locally isomorphic to the pure Fibonacci chain, and therefore 
our results imply that locally isomorphic chains, besides having the same energy 
spectrum,^^ also have the same Lyapunov exponent, and hence their eigenstates 
have the same nature as the ones of the pure Fibonacci chain. We are now in- 
vestigating the effects of random tiling on the Lyapunov exponent, of electron 
and phonon excitations, on other aperiodic chains, such as the Thue-Morse, the 
period-doubling, and binary non-Pisot sequences. Other types of disorder are 
also being considered on the Fibonacci chain, as well as on the other aperiodic 
chains mentioned, in order to understand the relevance/irrelevance of disorder 
on the Lyapunov exponent, and consequently on the localization properties of 
those systems. The results of this work will be reported elsewhere. 
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Figure captions 

FIG. 1. Inverse of the phonon wavefunction at the Fibonacci sites n = i^i+i, 
for energy E = 0.4, on a pure Fibonacci chain, and the associated Lyapunov 
exponent 7. 

FIG. 2. Lyapunov exponent 7 for: a) electronic "diagonal" (e^ = — £/3 = s-y = 

1) , b) electronic "off-diagonal" {Va = 1, Vb = 2), and c) phonon (Va = 1, Vb = 

2) excitations on a pure Fibonacci chain. 

FIG. 3. Self-similar structure of 7: a) electronic "diagonal" (e^ = — £/3 = s-y = 

1) , b) electronic "off-diagonal" (Va = 1, Vb = 2), and c) phonon (Va = 1, Vb = 

2) excitations, to compare with FIG. 2. 

FIG. 4. Maximum 7 in a gap vs gap width AEg, for: a) electronic "diagonal", 
(A) (£„ = -Sfj = = 1, 5 = 0.62 ), (■) {sa = -ef3 = = 2, S = 0.47); b) 
electronic "off-diagonal", (•) (Va = 1, Vb = 2, 5 = 0.75), (♦) (Va = 3, Vs = 1, 

S = 0.53), excitations on a pure Fibonacci chain. 

FIG. 5. Maximum 7 in a gap vs gap width AEg, for phonon excitations (Va = 1, 
Vb = 2), on a pure Fibonacci chain. 

FIG. 6. Power-law exponent 6, of 7max vs AEg, for phonon excitations: a) 
Va = 1, Vb = 2, b) Va = 2, Vb = 1, on a pure Fibonacci chain. 

FIG. 7. Lyapunov exponent 7 vs probability of disorder p, for phonon excita- 
tions, with energies: (•) E = 1.2, (♦) E = 2.3, (a) £ = 3.4, on random periodic 
chains. 
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FIG. 8. Lyapunov exponent 7 vs probability of disorder p, for various energies 
E, of: a) electronic "diagonal", (•) E = -1.9, (♦) E = 0.15, (a) E = 1.1; b) 
electronic "off-diagonal", (•) E = 0.5, (♦) E = 1.5, (a) E = 2.05, c) phonon, 
(•) E = 1.4, (♦) E = 3.1, (a) E = 5.29, excitations on random Fibonacci 
chains. 
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